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Abstract 

The distance energy of a graph G is a recently developed energy- type invariant, defined as 
the sum of absolute values of the eigenvalues of the distance matrix of G. There was a vast 
research for the pairs and families of non-cospectral graphs having equal distance energy, 
and most of these constructions were based on the join of graphs. A graph is called circulant 
if it is Cayley graph on the circulant group, i.e. its adjacency matrix is circulant. A graph 
is called integral if all eigenvalues of its adjacency matrix are integers. Integral circulant 
graphs play an important role in modeling quantum spin networks supporting the perfect 
state transfer. In this paper, we characterize the distance spectra of integral circulant graphs 
and prove that these graphs have integral eigenvalues of distance matrix D. Furthermore, we 
calculate the distance spectra and distance energy of unitary Cayley graphs. In conclusion, 
we present two families of pairs (Gi,G2) of integral circulant graphs with equal distance 
energy - in the first family Gi is subgraph of G2 , while in the second family the diameter of 
both graphs is three. 

Keywords: distance matrix; distance energy; unitary Cayley graph; integral circulant 
graph. 
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1 Introduction 

Let G be a simple undirected graph with n vertices. The vertices of G are labeled as 0, 1, ... , n—1. 
The distance between the vertices i and j is the length of a shortest path between them, and is 
denoted by d{i,j). The diameter of G, denoted by diam{G), is the maximum distance between 
any pair of vertices of G. 

Let A be the adjacency matrix of G, and Ai, A2, . . . , A„ be the eigenvalues of the graph G. 
The energy of G is defined as the sum of absolute values of its eigenvalues [1] , 

n 

E{G) = ^\Xi\. 
1=1 

The energy is a graph parameter stemming from the Hiickel molecular orbital approximation 
for the total 7r-electron energy (for recent survey on molecular graph energy see [2j and [3J). 

The distance matrix of a graph G is the square matrix D{G) = [d{i, The eigenvalues 

of the distance matrix D{G), labeled as /ii ^ /.i2 ^ . . . ^ /^n, are said to be the distance or 
D-eigenvalues of G and to form the distance or L>-spectrum of G [1]. The sum of distance 
eigenvalues is zero, ^27=1 A** = 0- The characteristic polynomial and the eigenvalues of the 
distance matrix of a graph were considered in [5]-[llj. 



Definition 1.1 The distance energy DE[G) of a graph G is the sum of absolute values of the 
eigenvalues of the distance matrix of G. 

Distance energy is a useful molecular descriptor in QSPR modeling, as demonstrated by 
Consonni and Todeschini in [12]. To avoid trivial cases, we say that the graphs G and H of the 
same order are D-equienergetic if DE{G) = DE{H), while they have distinct spectra of distance 
matrices. 

Our motivation for this research came from various constructions of non-cospectral, equiener- 
getic graphs with equal number of vertices [13j-|16]. Indulal et al. [17] constructed pairs of 
D-equienergetic graphs on n vertices for n = 1 (mod 3) and for n = (mod 6). Ramane et al. 
[18] proved that if Gi and G2 are r-regular graphs on n vertices and diam{Gi) ^ 2 , i = 1,2, 
then DE{L^{Gi)) = DE{L''{G2)) for k ^ I, where L^{G) is the fc-th iterated line graph of G. 
In [12] the authors obtain the eigenvalues of the distance matrix of the join of two graphs whose 
diameter is less than or equal to two, and construct pairs of non D-cospectral, D-equienergetic 
graphs on n vertices for all n ^ 9. Stevanovic and Indulal [20J further generalized this result and 
described the distance spectrum and energy of the join-based compositions of regular graphs in 
terms of their adjacency spectrum. Those results are used to show that there exist a number of 
families of sets of non-cospectral graphs with equal distance energy, such that for any n £ N, 
each family contains a set with at least n graphs. All these constructions from the literature are 
based on graph products and most of presented graphs have diameter two. 

A graph is called circulant if it is Cayley graph on the circulant group, i.e. its adjacency 
matrix is circulant. A graph is called integral if all eigenvalues of its adjacency matrix are 
integers. Integral graphs are extensively studied in the literature and there was a vast research 
for specific classes of graphs with integral spectrum |21] . 

Integral circulant graphs were imposed as potential candidates for modeling quantum spin 
networks with periodic dynamics. For the certain quantum spin system, the necessary condition 
for the existence of perfect state transfer in qubit networks is the periodicity of the system 
dynamics (see [22]). Relevant results on this topic were given in [23], where it was shown that 
a quantum network topology based on the regular graph with at least four distinct eigenvalues 
is periodic if and only if it is integral. Various properties of integral circulant graphs were 
investigated in [2l]-[28]. 

Integral circulant graphs arise as a generalization of unitary Caley graphs, recently studied 
by Klotz and Sander [25]. Let D be a set of positive, proper divisors of the integer n > 1. Define 
the graph ICGn{D) to have vertex set Z„ = {0, 1, . . . , n — 1} and edge set 

E{ICGn{D)) = {{a, 6} I a, 6 e Z„, gcd(a -b,n)(^D}. 

In this paper our intention is to move a step towards in the investigation of the graph 
theoretical properties of integral circulant graphs that are important parameters of quantum 
networks. Namely, we deal with the distance matrix of integral circulant graph ICGn{D) and 
distance spectra. The integral circulant graph /CGio(l) is shown on Figure 1, together with its 
distance matrix. 

The paper is organized as follows. In Section 2 we present some preliminary results on 
integral circulant graphs, while in Section 3 we prove that distance matrix of ICG graphs have 
integral spectra, and characterize the vertices at distance k from the starting vertex 0. Our main 
result is the description of the distance spectrum of ICGn{D). In Section 4 we calculate the 
distance energy, the distance spectral radius and the Wiener index of unitary Cayley graphs. 
Note that the distance energy of ICGn(l), where n is even with odd prime divisor, is not 
fully resolved since the sign of expression 2{ip{n) — 1) — ^ can vary. Finally, in Section 5 we 
present two families of non-cospectral distance equienergetic graphs (/C G^p ( 1 ) , IG G^p ( 1 , p) ) and 
{IGG2pq{l,p), IGG2pq{l,q)) for arbitrary primes p > q > 3. These constructions are not based 
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Figure 1: The integral circulant graph ICGio{l) and its distance matrix. 

on the graph products, and in the first pair we have that ICGsp{l) is subgraph of ICG-ip{l,p), 
while in the second pair the diameter of both graphs is three. 

2 Preliminaries 

Let us recall that for a positive integer n and subset 5 C {0, 1, 2, . . . , re — 1}, the circulant graph 
G(re, S) is the graph with re vertices, labeled with integers modulo re, such that each vertex i is 
adjacent to |5| other vertices {i + s (mod re) | s € S}. The set S is called a symbol of G(re, S). 
So [23] has characterized integral circulant graphs. Let 

Gn{d) = {k I gcd{k,n) = d, 1 < k < n} 

be the set of all positive integers less than re having the same greatest common divisor d with 
re. Let Dn be the set of positive divisors d of re, with d ^ ^. 

Theorem 2.1 A circulant graph G{n,S) is integral if and only if 

S=(j Gn{d) 
d&D 

for some set of divisors D C Dn- 

It follows that the degree of IGGn{D) is equal to J2deD ^i''^/^)^ where <^(re) denotes the 
Euler phi function. Authors in [28] proved that an integral circulant graph ICGn{D) = 
IGGn{di,d2, . . . , dk) is connected if and only if gcd(di, (i2, • • • , dk) = 1. 

Let D be an arbitrary set of divisors {di, d2, ■ ■ ■ , d^} of re. We establish the distance matrix 
Dn{D) of integral circulant graph ICGn{D) with respect to the natural order of the vertices 
0, 1, . . . , re — 1. The entries oq, ai, . . . , a„_i of the first row of Dn{D) generate the entries of the 
other rows by a cyclic shift. 



DniD) 



( ao a\ ... On-i \ 

On-l OO • • • 0"a-2 



\ tJi 02 ... ao / 

For more details on circulant matrices see [29j . There is an explicit formula for the eigenvalues 
//.r,0^r^re— 1, ofa circulant matrix such as D„- Define the polynomial -Pn(-2) by the entries 



of the first row of D^., 

n-l 



j=0 
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The eigenvalues of -D„ are given by 

n-l 

l^r = Pniu^l = X] • ^ r ^ n - 1. (1) 

j=0 

Ramanujan's sum, usually denoted c(r, n), is a function of two positive integer variables r 
and n defined by the formula 



a—l a— 1 

gcd(a,n) — 1 gcd{a,n) — 1 



where Un denotes a complex primitive n-th root of unity. These sums take only integral values, 



c{r, n) = fi 



5cd(r, n) 



^ I gcd(r,n) 



where n denotes the Mobious function. In the next section, we will use the well-known summa- 
tion [SU] 

n-l 



if r \ n 
n if r \ n 



i=0 

In [25] it was proven that gcd-graphs (the same term as integral circulant graphs ICGn{D)) 
have integral spectrum, 

Afc = c (A;, ^) , (2) 

deD 

3 Distance spectra of ICG 

Let Np(i) be the set of vertices that are at distance p from the starting vertex i. We will prove 
that 

Sp 

iVp(0) = |jGn(d?^ 
i=l 

for some divisor set of n 

For p = 1, this fits the definition of integral circulant graph ICGn{D) for divisor set D^^'^ = 
D. Let a be an arbitrary vertex at the distance p from the starting vertex 0. Assume that 
gcd(a, n) = d. If follows that there exist vertices = vo,vi, . . . ,Vp-i,Vp = a, such that for all 
^ r ^ p — 1, the vertices Vp and t'p+i are adjacent. In other words, 

gcd(t'r. — Vr+i,n) = di^., for r = 0, 1, . . . ,p — 1, 1 ^ ir ^ k. 

Let b be an arbitrary vertex (different than a) such that gcd(6, n) = d. This means that 
there are integers a' and h' relatively prime with n, such that a = a' ■ d and b = b' ■ d. Since 
gcd(a',n) = 1, from Bezout's identity we have 

a ■ a* = 1 (mod n), 

where a* is the inverse of a' modulo n. Now, consider the following vertices = uq, ui, . . . , Up-i,Up 
b, defined as 

Ur = Vr • a*b'. 
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For r = we have uq = 0, and for r = p we have that Up = a ■ a*b' = d ■ {a' ■ a*) ■ b' = d ■ b' = b. 
For r = 0, 1, ... ,1? — 1 holds 

gcd{ur - Ur+i,n) = gcd{{vr - fr+i) ' a*b',n) = gcd{vr — Vr+i,n) = di^. 

This proves that the distance from to the vertex b is less than or equal to k, d{0,b) ^ d{0,a). 
Similarly, we prove that d{0, a) ^ d{0, b) and finally, we get d(0, a) = d{0, b) = k. 

It is easy to see that the set of all divisors is the union of all sets D(p\ p = 1, 2, . . . , diam{G). 
Using the relation ([1]), we get the distance spectra of ICGn{D): 

This proves the main result of this section. 

Theorem 3.1 Integral circulant graph ICGn{D), where D is an arbitrary set of divisors of n, 
has integral distance spectra. 

The Wiener index of G is the sum of distances between all pairs of vertices, 

W{G)= dia,b). 

a,beV 

The Wiener index is considered as one of the most used topological indices with high correlation 
with many physical and chemical properties of molecular compounds (for recent results and 
applications of Wiener index see [31] ) . 

In [TDj it is proven that //q ^ with equality if and only if all row sums are equal. 

Since the distance matrix of IGGn{D) is also a circulant matrix, it follows that 

wia) = 

4 Distance energy of UCG 

The Unitary Cayley graph IGGn{^) = Xn has the vertex set V{Xn) = Zn and the edge set 

E{Xn) = {(a, 6) : a, 6 € Z„, gcd(a -b,n) = 1}. 

The graph X„ is regular of degree ip{n). Unitary Cayley graphs are highly symmetric and have 
some remarkable properties connecting graph theory and number theory. Fuchs [32j showed 
that the maximal length of an induced cycle in X„ is 2'^ + 2, where k is the number of different 
prime divisors of n. Klotz and Sander determined the diameter, clique number, chromatic 
number and eigenvalues of unitary Cayley graphs. 

In next four subsections, we will calculate the distance spectra and distance energy of Xn- 

4.1 n is a prime number 

In this case, Xn is a complete graph Kn with diameter 1. The spectra of X„ is {p—1, — 1,— 1,...,— 1} 
and 

DE{Xn) = 2(n - 1). 
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4.2 n is the power of 2 

If n = 2^^' for k > 1, then X„ is a complete bipartite graph with the vertex partition 

ViXn) = {0, 2, . . . , n - 2} U {1, 3, . . . , n - 1}, 
and has diameter 2. The adjacency spectrum of is — ^, 0, 0, . . . , 0}, and consequently the 

in 
2 

' " ' n 



distance spectrum of X„ is - 2, f - 2, -2, -2,..., -2}. Since 93(2'=) = 2^="^ it follows that 



DE{Xn) 



3n 

2 

2 



+ 



2 
2 



+ 2(n - 2) = 4(n - 2). 



4.3 n is odd composite number 

We need the following result from 



Theorem 4.1 The number of common neighbors of distinct vertices a and b in the unitary 
Cayley graph Xn is given by Fn{a — b), where Fn{s) is defined as 



Since 2 does not divide n, according to Theorem 14. II all factors in the expansion of Fn{a — b) 
are greater than zero. Therefore, in this case there is a common neighbor of every pair of distinct 
vertices, which implies that the diameter is 2 and 

Dn{Xn) = 2(J„ — In) — j4„(X„), 

where Jn is the matrix of ones and In is the identity matrix. The eigenvalues of Dn are given 
by 

n-l 



//, = 2^<-2-c(r,n). 

1=0 

Finally, using tliG relation the spectra of Dn i^ 

{2(n - 1) - ip{n), -2 - c(l, n), -2 - c(2, n), . . . , -2 - c(n - 1, n)} 

Therefore, the distance spectral radius of Xn is 2(n — 1) — ^{n). We can calculate the distance 
energy of Xn using similar technique as in [33] ■ The nullity of a graph G, denoted as r]{G), is 
the multiplicity of zero as the eigenvalue in adjacency spectra. 

Lemma 4.2 The nullity of Xn is n — m, where m = piP2 ■ ■ ■ ■ • Pk is the maximal square-free 
divisor of n. 

Consider the following sum 

n n 

5" = J] I - c(i, n)-2\ = Y^ |c(i, n) + 2|. 
1=1 1=1 

We already know the nullity of Xn, and therefore we will sum only the non-zero eigenvalues 
from the spectra of Xn- Divide the sum S in two parts: when ^^^"^ is a square- free number 
with an even number of divisors and when gg^"„ is a square-free number with an odd number 
of divisors. The number of even subsets of {pi,P2, ■ ■ ■ ,Pk} equal to the number of odd subsets 
of {pl,^'2,•••,Pfc}• 
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In the first case, we have 

y:f-/^ + 2)=^(n)-2^-^ + 2 (3) 

i&Si \^^Scdit,n)) ) i|m, /i(/)=l 

Let I be a square-free number that divides m with an even number of prime factors. The 
number of solutions of the equation ^^^"^ = I is equal to </'(Z). For all ^ z < n that satisfy 
n = I • gcd(i, n) we have 



gcd(j,n) ' 



After taking the summation for all / | m with = 1 we derive the identity ([3]). Analogously, 
in the second case we have 

y:f-7^-2)=^(..).2^-'-2 Y, fo- 

Since Euler function ip{n) is multiplicative, after adding the above sums and eigenvalue —2 
we get 

S = 2 ■ {n - m) + ^{n) ■ 2'' + 2J2 

l\m 

k k 
= 2n-2m + ip{n) -2^ + 2 JJ(1 - ip{pi)) = 2n-2m + Lp{n) -2^+2 JJ(2 - pi). 



i=l i=l 

Finally, the distance energy of Xn equals 

DE{Xn) = 5- |(/9(n) + 2| + |2n-(^(n) -2| 

= 2 ( 2n + v3(n)(2''~^ -1) -m 



k \ 

2+n(2-p.)) 

i=i J 



4.4 n is even with odd prime divisor 

As in the previous cases there are ^{n) ones in the first row of matrix D„. The vertex is 
not adjacent to even vertices 2, 4, . . . , n — 2, but any two even vertices have a common neighbor 
by Theorem 14.11 Since all neighbors of an odd vertex are even vertices, we conclude that the 
number of vertices at distance 2 is exactly § — 1- Similarly, any two vertices a and h of X„, 
which are both odd have a common neighbor. This proves that the diameter of Xn is 3 and 
there are exactly ^ — (p{n) vertices at distance 3 from the starting vertex 0. 
It follows that the polynomial Pn[z) equals 

1 ——1 

n— 1 2 n 

P„(z) = 3^z' - ^^2^-2-2 ^ z\ 

i=0 i=0 »=l 

gcd(i,ri) — 1 

The eigenvalues of -D„ are given by 

1 -—1 

n—l 2 

/x, = 35]<-^.;^-2-2.c(r,n). 

1=0 i=Q 
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For r = 0, the distance spectral radius is equal to 



/io = 3n - ^ - 2 - 2^{n) = ^ - 2{ip{n) + 1). 



For r = ^, it follows that c{n/2,n) = —ip[n) and 



n 



,n 



n 



Mn/2 = 0---2-2.c(-,n) = 2(<^(n)-l) ^ 



For r 7^ and r 7^ S, the value of the sum 



-1 



i=0 



is zero. Finally, the spectra of Dn consists of ^ — 2{(p{n) + 1), 2{ip{n) — 1) — ^ and 

-2 -2 -0(1, n), -2 -2 -0(2,71), . . . ,-2-2 -0(^/2-1,71), -2-2 -0(71/2 + 1,71), . . . - 2-2-c(7i-l,7i). 

The sum S' = Y^I^q |c(i,7i) + 1| can be computed analogously, and since the smallest prime 
of 71 is = 2 we have 



S' = n-m + ip{n) ■ 2^ + JJ(2 - pi) = n - m + (p{n) - 2^ 

i=l 

The distance energy of X„ is equal to 
DE{Xn) = 25' -|2 + 2-c(0,7i)|-|2 + 2-c(7i/2,7i)|+ ^-2(^(77) + 1) 

/ 571 



+ 



2{^{n) 



n 
2 



271 - 2771 + ip{n) ■ 2''+^ - (2 + 299(71)) - (2(^(n) " 2) + ( y - 2((/?(n) + i; 



+ 



71 



2{^{n) - 1) - 2 



The value of |2((^(7i) — 1) — §1 cannot be resolved, since it takes both positive, zero and 
negative values (consider for example 71 = 6, 71 = 10 and 71 = 12). By using the representation 
71 = 2'^ -771 with 771 being odd, it follows '^{n) = Lp{2^)(f{m) = 2^~^ip{m). Therefore, the equation 
2{ip{n) — 1) = ^ is equivalent with 

2v3(7l) = 71+ 1. (4) 

The equation (4) is related to the still open conjecture of Lehmer [3l], and some obvious solutions 
involve prime Fermat numbers, n = Ff^ — 2 = 2 — 1 for A: = 0, 1, 2, 3, 4, 5. 



5 Distance equienergetic graphs 

Let 71 = pq, for odd prime numbers p and q. According to Subsection 4.3, the distance energy 
of /CG„(1) equals 

DE{ICGn{l)) = 3pq + 2{p-l){q-l)-n + 2{2-p){2-q) 
= 6{p-l)iq-l). 

Consider the graph ICGn{^,p) ~ its diameter is two, and only the vertices satisfying gcd(a — 
b,n) = q are at distance 2. Using the formula ([2]) we get the distance eigenvalues of ICGn{^,p), 

fir = c{r,pq) + c(r, q) + 2 ■ c{r,p). 
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Based on the greatest common divisor of r and pq, it follows that 

1-1 + 2- (-1) ifpfr and q\r 



fir 



-ipip) - 1 + 2 • ifip) 
(p{pq) + (p{q) + 2 • ip{p) 



if p I r and qjr 
if p I r and q \ r 
if p I r and q \ r 



Finally, the spectra of ICGn{l,p) consists of —2 with multipUcity pq — q,p — 2 with multiplicity 
g — 1, and pq + p — 2 with multiplicity 1. Therefore, the distance energy equals 



DE{ICGn{l,p)) 



= q{p-l)-2 + {q 
= 4g(p-l). 



l).(p-2) + l-{pq+p-2) 



Only for g = 3 and n = Sp, we have the identity DE{ICGn{l)) = DE{ICGn{l,p)), which 
represents a family of distance equienergetic integral circulant graphs. Note that /CG3p(l) is a 
subgraph of ICG3p{l,p) for all primes p > 3. 

Let now n = 2pq, where p and q are different odd prime numbers. We will prove that graphs 
ICGn{l,p) and ICGn{l,q) have equal distance energy, 

DE{ICGn{l,p)) = DE{ICGn{l,q)) = 4.{3pq ~ p - q - 1). 

According to Subsection 4.4, the distance eigenvalues of bipartite ICGn{^,p) are given by 

fij. = c{r, 2pq) + c(r, 2q) + 2c{r,pq) + 2c{r,p) + 2c(r, q) + 3c(r, 2p) + 3c(r, 2). 

Based on the greatest common divisor of r and 2pq, it follows that 



-1+1+2-2-2+3-3 
1-1+2-2-2-3+3 

(p - 1) + 1 - 2(p - 1) + 2(p - 1) - 2 - 3(p - 1) - 3 

-(p - 1) - 1 - 2(p - 1) + 2(p - 1) - 2 + 3(p - 1) + 3 
(g - 1) - (g - 1) - 2(g - 1) - 2 + 2(g - 1) + 3 - 3 
-{q - 1) + (g - 1) - 2(g - 1) - 2 + 2(g - 1) - 3 + 3 
pq-2p-2 
3pq + 2p-2 



if 2 I r and p | r and q \ r 

if 2 I r and p \ r and g | r 

if 2 f r and p \ r and g f r 

if 2 I r and p \ r and q \ r 

if 2 I r and p f r and g | r 

if 2 I r and p \ r and g | r 

if 2 I r and p \ r and g | r 
if 2 I r and p \ r 



and 



3pq + 2p-2 

pq — 2p — 2 
-2p-2 
2p-2 
-2 



Since the multiplicity of distance eigenvalues 2p 

DE{ICGn{l,p)) = {3pq + 2p-2) + {pq-2p 
= \2pq - 4p - 4g - 4, 



for r = 

for r = pq 
for gcd(r, 2g) 
for gcd(r, g) = 
otherwise 

2 and 



1 



2p — 2 is g — 1, we finally get 

2) + (g - l){2p + 2) + (g - l){2p - 2) + 2{pq 



2q) 



which is a symmetric expression involving p and g. Note that the diameter of graphs ICG2pq{l,p) 
and ICG2pq{l,q) is equal to three. 

These constructions of infinite families of distance equienergetic graphs arc the first ones 
derived not using the product of graphs nor iterated line graphs. It would be of great interest 
to calculate the distance energy of other classes of integral circulant graphs, and explore new 
families of non-cospectral distance equienergetic integral graphs. 
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